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Простой алгоритм цифровой 
параметрической идентификации 


Рассматривается цифровой алгоритм определения коэффициента передачи по угловой скорости и обеих 
постоянных времени математической модели электропривода, когда его передаточную функцию по углу 
можно представить в виде последовательного соединения передаточных функций двух инерционных и одного 
интегрирующего звеньев. Описана структурная схема алгоритма, представлены результаты моделирования. 


Введение 


Наличие математической модели (ММ) объекта управления (ОУ), описывающей 
его как можно более точно, является важным условием разработки эффективных ал- 
горитмов управления современными системами автоматического управления (САУ). 
При этом весьма желательным и важным является знание изменения параметров мо- 
дели в зависимости от условий эксплуатации САУ (изменение температуры, влажности, 
присутствия агрессивной среды и т.п.). 

В терминах теории анализа электрических цепей задача идентификации достаточно 
исследована и решение ее в основном сводится к анализу переходной функции электри- 
ческой цепи при подаче на ее вход единичного скачка тестирующего напряжения [1]. 

На этом принципе построены известные процедуры идентификации [2], которые 
сводятся к анализу напряжения 2({) на выходе датчика углового положения, который 


жестко механически связан с выходным валом электропривода (ЭП), при воздействии 
на вход разомкнутого электропривода тестирующего напряжения х(). Однако в этом 
случае возможна идентификация только самой простейшей ММ электропривода. 

В [1] [3] предложены алгоритмы идентификации для определения коэффициента 
передачи по угловой скорости А и постоянных времени 7, Г, для более сложной ММ 


привода, когда его передаточную функцию по углу И’икр) можно представить в виде 
последовательного соединения передаточных функций двух инерционных и одного 
интегрирующего звеньев. Эти алгоритмы позволяют с достаточной для практического 
применения точностью определять коэффициент передачи по угловой скорости электро- 
привода и обе его постоянные времени. Однако для их реализации необходимо исполь- 
зовать достаточно сложные итерационные алгоритмы формирования специальной ап- 
проксимирующей функции переходной характеристики по скорости тестируемого ЭП. 

Целью данной работы является анализ более простого алгоритма цифровой пара- 
метрической идентификации, который позволяет с достаточной для практики точностью 
определять коэффициент передачи по угловой скорости А и постоянные времени Т/, Т> 
ММ электропривода, когда его передаточную функцию по углу можно представить в 
виде последовательного соединения передаточных функций двух инерционных и одного 
интегрирующего звеньев. 
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Структурная схема алгоритма представлена на рис. 1. 


О ИИА ИТ А, 


Убтт ет 


Рисунок 1 — Структурная схема эффективного алгоритма цифровой параметрической 
идентификации электропривода 

На рис. 1 принято, что х(Ё) — скачок напряжения с амплитудой А на входе ЭП; 

2(©) — отклик аналогового датчика углового положения (ДУП) с передаточной функцией 

И/луп(р) на х(0 (переходная характеристика (ПХ) ЭП по угловому положению); 2(иТ,) — 

отсчеты ПХ ЭП на выходе аналого-цифрового преобразователя (АЦП); у(Р) — отклик 

аналогового датчика угловой скорости (ДУС) с передаточной функцией И’, „с (1) на х(0) 


* 


(ПХ ЭП по скорости), у(иТ,) — отсчеты ПХ по скорости; 1, 1. — оценки первой и 


* 
второй постоянных времени ЭП; К -— оценка коэффициента передачи по угловой ско- 
* 
рости ЭП; ВСПК - блок вычисление оценки т, суммы постоянных времени т, = 1 +Т, 


‚ 
и оценки К коэффициента А; ВВУ — блок вычисления номера отсчета времени пере- 
ходного процесса и’, соответствующего времени установления переходного процесса, 
т.е. времени при котором скорость вращения двигателя ЭП достигает значения 0,95 от 
своего номинального значения ©), и номера отсчета и! времени переходного процесса, 


при котором скорость вращения двигателя ЭП достигает значения 0,906 „„; ВПВ2 — 
блок вычисления оценки второй постоянной времени Т»; ВПВ1 — блок вычисления 


оценки постоянной времени т. : 

Процедура идентификации в соответствии с предлагаемым алгоритмом сводится к 
тому, что на вход разомкнутого ЭП подается скачок напряжения, а в ЭВМ выполняются 
вычисления в соответствии с процедурой, представленной на рис. 1. 

В [1] получена формула зависимости 2({) , где для простоты выкладок полагалось, 
что передаточная функция датчика углового положения И’дул (р) = 1. Перейдя к дискрет- 
ной записи этой формулы, можно показать, что поведение 2(пТ, ) описывается прибли- 
женной формулой: 


2 2 


—Т пт Т пт 
2(пТ,) = О 1_ех °)+—2 ех уе ИГ. 1 
( 9) ном т = р( Т ) Г, ЕЙ р( Е ) ил 1 д ( ) 
где О, =К-А, и — номер отсчета входного сигнала, Г, — период дискретизации, для ко- 


> п 
торого в соответствии с теоремои Котельникова должно выполняться условие: Т д А 
0) 


с 


где @, — частота среза системы. 


В [! также получена формула для уравнения асимптоты Р({) функции 2(9. Перейдя 
к дискретной форме записи этого уравнения, получим: 


Р(иТ,) = Р(иТ,) |т,„=О „ИТ, — СТ, -Т,)} (2) 


пу — 
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Из последнего соотношения следует, что если Р(иТ,) определено и А известно, 


то могут быть найдены оценки 7, и’. 

Дополнительные исследования показали, что с методической ошибкой, не превы- 
шающей нескольких процентов, для определения т) и К” вместо уравнения асимптоты 
достаточно по результатам обработки 2(иТ,) находить уравнение касательной к функции 
2("Т,) в точке (Пт, 2(И т) )- 

В соответствии с вышесказанным в блоке ВСПК для нахождения оценок т) И 
К’ последовательно выполняются следующие действия: 

— по результатам обработки текущих значений 2(иТ,) на всем интервале наблюде- 


ния находят уравнение касательной к функции 2(пТ,) в точке (и„Т», 2(п„Г,)), кото- 


уст уст 


рое в общем виде имеет вид: Р(иТ,) = В\иТ, -С}; 
* * * В 
— находят искомые значения т, и К по формулам: А = р т, =С. 


В [3] также получена формула для описания зависимости у(!) , где для простоты вы- 
кладок полагалось, что передаточная функция датчика углового положения И’ус(р)=1. 
Перейдя к дискретной записи этой формулы, можно показать, что поведение у(иТ,) опи- 


сывается приближенным соотношением: 


Т, ПТ Т. пт 
ПТ) =О те о 2 ех Ех, Э 
У( и ном т, — р р( Т. ) Т, т р( г ) ( ) 


из которого следует, что, сохранив в памяти ЭВМ значения у(иТ,) в диапазоне ско- 


* 


ростей, лежащих от нуля до её установившегося значения, можно, кроме © 


* * 
делить оценки номеров отсчета времени П уст ип. 


опре- 


ном ? 


В соответствии с этим в блоке ВВУ по результатам экспериментально полу- 
ченных значений у(иТ,) последовательно находятся значения ©”, 0,9550; п 


ном ? уст 2 
0,900 


* * 


п, . Затем значения и„„ и и, передаются в блок ВПВ2. 


ном ? ‘уст 


В блоке ВПВ?2 производится вычисление оценки второй постоянной времени Т. — 


Процедура получения оценки базируется на том факте, что при существенном неравенстве 
значений постоянных времени 7! и 1> переходным процессом, вызванным составляющей 
с меньшей постоянной времени при приближении разгонной кривой к своему установив- 
шемуся значению, можно пренебречь. Тогда большую постоянную времени можно оце- 
нить, например, по разности времени Д! между временем разгона двигателя ЭП до 


* 


скорости 0,95 О. и временем его разгона до скорости 0,90 0 


ном * 


Действительно, как показала практика использования алгоритма, при выполне- 


п уст д Т, 


нии условия пГ, > = и Т2> 1,7 Т,в (4) составляющей 


1 
ехр( й ) суммы 
2 1 1 

в фигурных скобках (1) можно пренебречь и на основании этого записать систему из 
двух уравнений: 


п. То 


Т. 
у(пт.Ть) 2 =0 1 2 ехр( 


а“ 2 
т 


) = 0,902, (4) 


1 2 
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Т. И уст То 
У(пТЬ) трет ы ном — ехр( т ) = 0,95 ° эт ы (5) 


2 1 2 
В результате решения системы из двух уравнений получим формулу для опре- 
деления значения второй постоянной времени: 


Ти -п 
= и =1,443-Т,(п„-п). (6) 


и. 
1-0,95 
Отсюда следует, что процедура определения Т› в блоке ВПВ2 сводится к вычисле- 


ниям по формуле (6), где вместо п„„ и и, необходимо соответственно использовать их 


уст 


оценки и„ии,. 


уст 


В блоке ВПВ1 7” рассчитывается по формуле: 
* * * 
1=1, -Т,. (7) 
Для проверки работоспособности алгоритма в пакете программ «МаЧа6» проводи- 
лось математическое моделирование его работы для имитационной модели тестируемого 
электропривода. 
Для этого в имитационной модели ЭП задавались постоянные времени Ть, Т> и 
коэффициент передачи по скорости А, а также амплитуда А тестирующего напряжения. 
В процессе моделирования алгоритма проводились три эксперимента, где Т> 
было выбрано равным 0,5 с, а Т, задавалась равной соответственно 0,05 с, 0,2 с и 


0,3 с. Амплитуда А задавалась равной / В, К = 5 и ‚ а УИ = 0,02 с. 
[@ . 


Результаты моделирования представлены в табл. 1. 


Таблица 1 
Эксперимент 1 Эксперимент 2 Эксперимент 3 
Г, =0,050 |Т,=0,500 | Т,=0,200 | Г, =0,500 |Т,=0,300 | Г, =0,500 
Т'=0,0425 | Т,=0,501 | Т/=0,195 | Т, =0,503 | Т,=0,289 | Т,=0,521 
71 =5% 72=2.0% [71=4,5% [|у2=6% 71 = 3,1% 72 =4,2% 


В табл. 1 принято, что у; и > — относительная погрешность измерения, соответст- 
венно, первой и второй постоянных времени. 

Результаты моделирования показали, что погрешность измерений постоянных вре- 
мени алгоритма достаточно высокая и не превышает 15% во всем диапазоне проведенных 
исследований, а погрешность определения коэффициента передачи по скорости привода 
не превысила 2%. 


Выводы 


Результаты проведенных исследований и имитационного моделирования позволяют 
сделать обоснованный вывод о том, что предложенный цифровой алгоритм параметри- 
ческой идентификации электропривода достаточно прост и эффективен. Его использование 
позволяет с достаточной для инженерного применения точностью определять коэффициент 
передачи по угловой скорости и обе постоянные времени электропривода для случая, 
когда его математическая модель может быть достоверно представлена в виде последова- 
тельного соединения передаточных функций двух инерционных и одного интегрирую- 
щего звеньев. 
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Простий алгоритм цифрово! параметрично! {дентифкацИ 

Розглядаеться цифровий алгоритм визначення коефицента передач! за кутовою швидюстю та обох 
сталих часу математично! модел! електропривода, коли його передавальну функщю за кутом можливо 
представити у вигляд! послдовного з’еднання передавально! функций двох 1нершиних та одне 1нтегруючот 
ланок. Описана структурна схема алгоритму, представлен! результати моделювання. 
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